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It is shown that a two-qubit phase gate and SWAP operation between ground states of cold trapped 
ions can be realised in one step by simultaneously applying two laser fields. Cooling during gate 
operations is possible without perturbing the computation and the scheme does not require a second 
ion species for sympathetic cooling. On the contrary, the cooling lasers even stabilise the desired time 
evolution of the system. This affords gate operation times of nearly the same order of magnitude 
as the inverse coupling constant of the ions to a common vibrational mode. 
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For many years now ion trapping technology has been 
one of the standard techniques for investigating quantum 
phenomena with single particles. Paul traps or more 
complex structures have been used to create chains of 
several hundred ions [1]. Many schemes for implementing 
quantum gates using trapped ions have been proposed. 
Some of them require cooling into the vibrational ground 
state [2-5] . The feasibility of this has already been exper- 
imentally demonstrated [6]. Other schemes can even be 
implemented with "hot" ions [7] and have been applied 
to entangle up to four ions [8] and to observe violation of 
Bell's inequality with single atoms [9]. 

This paper proposes an alternative scheme for realising 
a phase gate and SWAP operation with cold trapped ions 
in one step by simultaneously applying two laser fields. 
Compared with other schemes [2,4], the experimental ef- 
fort for quantum computing can thus be greatly reduced. 
Ideally, the system remains during the whole gate opera- 
tion in the ground state of a common vibrational mode. 
Thus cooling of this mode does not perturb the computa- 
tion. On the contrary, it is shown that cooling can even 
improve the gate fidelity significantly. This affords gate 
durations of nearly the same order of magnitude as the 
inverse coupling constant of the ions to the common vi- 
brational mode and there is no need to include a different 
ion species for sympathetic cooling [10] in the trap. In 
addition, the scheme only requires good control over one 
Rabi frequency. The size of the coupling strength of the 
ions to the vibrational mode does not enter the effective 
time evolution of the qubits. 

Each qubit is obtained from the atomic ground states 
|0) and 1 1) of one ion, while the common vibrational mode 
is cooled down to its ground state |0 V ib)- To establish 
coupling between qubits, we use a metastable state |2) 
and a strong laser field with Rabi frequency Q2 detuned 
by the phonon frequency v, as shown in Figure 1. The 
phase gate that adds a minus sign to the amplitude of 
the qubit state 1 01) then requires in addition only a weak 
laser pulse with the Rabi frequency fl coupling resonantly 



to the 0-2 transition of ion 1. Let us denote the creation 
operator of a single phonon in the mode as and intro- 
duce the coupling constant of the ions to the common vi- 
brational mode as g 2 = 5^2- (Here r/ is the Lamb-Dicke 
parameter characterising the steepness of the trap.) The 
Hamiltonian of the system within the dipole and the ro- 
tating wave approximation and in the interaction picture 
with respect to the free Hamiltonian is then given by 
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H = J2 i%2 |l)i(2| ftt + i m |0> x <2| + h.c. (1) 

i=l 



Here the Lamb-Dicke regime and the condition g\ <C v 2 
have been assumed, as in [2]. 
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FIG. 1. Level scheme for implementation of a phase gate. 
A strong laser field with detuning v establishes coupling be- 
tween the ions via a common vibrational mode. In addition, 
a laser pulse individually addressing the 0-2 transition of ion 
1 is required. 

Let us first assume that the coupling constant g 2 is a 
few orders of magnitude larger than the Rabi frequency 
Q. Then there are two different time scales in the system 
and the time evolution can be calculated to very good 
approximation by adiabatic elimination. To do so, the 
amplitude of the state with n phonons in the vibrational 
mode and the ions in is denoted as c n ij. Only the co- 
efficients of the qubit states, coocb cooi 7 C010 and con, and 
of the entangled state \0a) with |a> = (|12) — |21>)/v^ 
change slowly in time. Their time evolution is given by 

cooo = — 5 ^ c ° 20 ' 

COOl = T^7| ^ (cQa — Cos) , 
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CQa = " Cool 

coin = con = . (2) 

Setting the derivatives of all other coefficients equal to 
zero yields C020 = c 0s = 0. For fl <C 32, the time evo- 
lution of the system (2) can thus be summarised in the 
effective Hamiltonian 

H cS = -^M}[|001)(0a|+h.c] . (3) 

If the duration of the laser pulse equals T — 2v / 2~7r/ r 2, 
then the resulting evolution is the desired phase gate. 
As the effective Hamiltonian H c g and the gate operation 
time T are independent of the coupling constant c/2, the 
proposed scheme is widely protected against fluctuations 
of this system parameter. 

Deviations from the time evolution (3) arise because 
population accumulates unintentionally in the states 
1 110) and I 111). In first order in f2/#2 one has 

in m 

cno — — 7; — c ooo , cm — — - — cooi (4) 
and the fidelity of the proposed phase gate equals 



f(t, m = 1 - 




Figure 2 compares this fidelity with the exact solution re- 
sulting from numerical integration of the time evolution 
(1). Good agreement is only found for -C 32- For some- 
what larger Rabi frequencies, in general fidelities worse 
than the result predicted by (5) are obtained due to non- 
adiabaticity. Thus the fidelity is above 99 % for all initial 
states if il < 0.1 32 is chosen. This corresponds to gate 
operation times T > 90/<72- In the following we aim at 
enlarging the parameter regime for which the fidelity is 
at least 99% [11]. 
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FIG. 2. The fidelity of a single phase gate as a function 
of f2/(?2 f° r the initial qubit states 1 00) (a,b) and 1 01) (c,d). 
The curves (a,c) result from an exact solution of the time 
evolution (1) while the curves (b,d) result from (5). For the 
states 1 10) and |11) optimal fidelities F = 1 are obtained. 

Dominating error source in the scheme is heating. To 
avoid this, the gate should be performed fast. Therefore 



we assume in the following that f2 is of nearly the same 
size as the coupling constant <?2- As can be seen from 
above, increasing f2 leads to the population of unwanted 
states. To reduce the error rate of the scheme one could 
therefore measure the population of states with n > at 
the end of each gate. Under the condition that no popu- 
lation is found in these states, the system gets projected 
back onto the subspace with n = and the fidelity of the 
prepared state increases. In case of the detection of an 
error, the gate failed and the whole computation has to 
be repeated. 

To realise such an error detection measurement one 
could use a laser field that couples the atomic ground 
state |1) with detuning v to an auxiliary state |3) and 
another laser that excites the atomic 1-2 transition with 
the same detuning. To assure that populating level 3 
leads to the emission of many photons, an even stronger 
laser should couple |3) to a rapidly decaying fourth level 
which decays into |3) with a high spontaneous decay rate. 
Gate failure leads thus to an effect which is known as a 
"macroscopic light period" [12], i.e. fast Rabi oscillations 
between level 3 and 4 accompanied by photon emission 
at a high rate (for details see [13]). This can easily be 
detected and the whole computation can be repeated, if 
necessary. 
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FIG. 3. Level scheme of a dissipation-assisted phase gate 
with cold trapped ions. In addition to the basic setup shown 
in Figure 1, a strong laser field couples the state |1) with 
detuning v to an auxiliary atomic state |3). Populating this 
level leads to the emission of photons with a rate IV 

However, such an error detection measurement would 
take much longer than the inverse of the coupling con- 
stants of the ions to the vibrational mode and much 
longer than T . Hence, error detection does not help to 
decrease the gate operation time for a given minimum fi- 
delity F. To shorten the gate duration without reducing 
the fidelity more than predicted by (5), we propose that 
the desired time evolution (3) be stabilised during the 
gate performance by using dissipation. This is achieved 
by continuously applying the laser fields, proposed for 
the implementation of error detection, as shown in Fig- 
ure 3. For simplicity, the strong laser coupling to the 
fourth level and spontaneous decay from this level have 
been combined into a single decay rate assigned to the 
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metastable state |3). Let us denote this spontaneous de- 
cay rate as r 3 , while g 3 = ^77^3 is the coupling strength 
of the atomic 1-3 transition to the vibrational mode. In 
the following, g 3 <~ g 2 and 10 g 2 < T 3 < 100 g 2 is as- 
sumed. 

More general, any process that indicates whether the 
phonon mode is excited or not can serve as an error detec- 
tion measurement. This applies to ground state cooling 
[14] because populating the vibrational mode leads to the 
emission of photons at a high rate while no emission takes 
place if the ions are in the vibrational ground state. Thus 
level 3 and the additional strong laser field in Figure 3 
could as well be replaced by the cooling laser setup. In- 
deed, continuous cooling can improve the fidelity of the 
performed gate operation. However, for simplicity, the 
continuous read out of the phonon mode is in the follow- 
ing modeled as shown in Figure 3. 

Basic mechanism of the improved scheme is that ob- 
serving for emitted photons implements a (conditional) 
no-photon time evolution, thus resulting in continuous 
damping of the population in unwanted states. As long 
as the amplitudes C020 and co s are negligible, the time 
evolution of the other states resembles the desired phase 
gate, as can be seen from (2). In addition, we show that 
the population that now accumulates in the states 1 1 10) 
and I 111) is about the same as predicted for an adiabatic 
process and the fidelity of the phase gate coincides with 
F(T, \tp)) given in (5) to a very good approximation. The 
price one has to pay for this improvement of the precision 
of the gate is that photon emission might occur with a 
small probability. Then the gate operation would have 
failed. 

To describe the time evolution of the system under the 
condition of no photon emission, we use in the following 
the Schrodingcr equation with the conditional Hamilto- 
nian i? CO nd- As predicted by the quantum jump approach 
[15], the norm of a vector developing with this non- 
Hermitian Hamiltonian decreases in general with time 
and 

Po(T,|^)) = ||C/ cond (T,0)|^)|| 2 (6) 

is the probability of no photon emission in (0, T) if |i/>) is 
the initial state of the system. For the level configuration 
shown in Figure 3 the conditional Hamiltonian equals 

2 3 

ffcond - E I^O'I b 1 + k m |0>i<2| + h.c. 

i=l j=2 

-^r 3 |3M3|. (7) 

Because of the different time scales of the scheme, the 
no-photon time evolution of the system can again be cal- 
culated by adiabatic elimination. This yields the same ef- 
fective Hamiltonian as in (3). In first order in S1/.92, pop- 
ulation accumulates unintentionally in the states |020), 
1 110), |030), |0s) |111), |013) and |031) and it is 



, > ifi V^gi l 53 93 \ 

(COs,Ciii,Coi3,C 03 l) — -- — I —■ -, — , — I Cool , 

2# 2 v 92^3 2 r 3 r 3 / 

I \ in ( S3 1 .93 \ , Q H 

(C020, ciio, c 03 o) = - — [^r, r~~ 3 ) Coo ° • (8) 

To optimise the fidelity, T 3 should be much larger than g 3 
so that all coefficients proportional g 3 /r 3 become negligi- 
ble. In this case, F(T, becomes the fidelity calculated 
in (5) to a very good approximation. 

That this is indeed the case is shown in Figure 4 which 
results from a numerical solution of the time evolution 
(7). As expected, the dissipation channel continuously 
introduced in the system stabilises the desired time evo- 
lution and corrects for errors resulting from the non- 
adiabaticity of the scheme if becomes of about the 
same order of magnitude as g 2 . The gate fidelity (5) ap- 
plies now to a much wider parameter regime. Fidelities 
above 99% are obtained if ft < 0.18 g 2 is chosen. How- 
ever, for too big spontaneous decay rates T 3 , the damping 
of unwanted amplitudes becomes ineffective which is why 
T 3 < 100 c/2 has been assumed. 
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FIG. 4. Fidelity of a single phase gate under the condition 
of no photon emission as a function of Q./g2 for the initial 
qubit states \ip) = |00) (a,b) and |V>) = |01) (c,d). Deviations 
from the time evolution with 53 = Ta — are continuously 
damped away with g A = 52 and 1^3 = 20 32 (a,c) and the 
fidelity is close to the theoretically predicted fidelity (5) as- 
suming adiabaticity (b,d). 

Using the coefficients co s and C020 given in (8) and the 
differential equations that govern the time evolution of 
the qubit states and the entangled state |0a) (the same 
as (2)), the unnormalised state of the system at the end 
of the gate operation under the condition of no photon 
emission can be calculated up to first order in fl/g 2 . Its 
norm squared equals the gate success rate (6) and 

p (r,|^)) = i-2V2^M[| C000 |2 + i| C001 | 2 ] _ (9) 

92 L 3 

This is in good agreement with the numerical results 
shown in Figure 5. For example, for £1 < 0.1 g 2 and 
T 3 = 20 g 2 one has Po(T,ip) > 95%, independent of 
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the initial state of the system. The probability for pho- 
ton emission during the gate operation is of the order of 
S1/.92 and for ^3/(5^3) <C 1 close to unity. Gate fail- 
ure might be a bit more likely than for quantum error 
detection (see Figure 5). However, no additional time 
is required which would increase the sensitivity of the 
scheme with respect to heating. 
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FIG. 5. Success rate of a single phase gate as a function 
of 0/<?2 for the initial qubit states \ip) — 1 00) (a,b) and 
\ip) = 1 01) (c,d) for the case of continuous monitoring of the 
system (a,c) and the case of only a single error detection 
measurement at the end of the gate (b,d). 

Another two-qubit gate that can easily be realised with 
the same experimental setup is SWAP operation. This 
gate exchanges the states of two qubits without the cor- 
responding ions having to change their places physically. 
Compared with the above phase gate, its implementation 
does not require individual laser addressing. To realise 
SWAP operation, the laser field with Rabi frequency Q, 
should address both ions simultaneously. To improve the 
gate fidelity, the same ideas as described above can be 
used. By analogy with (7), the conditional Hamiltonian 
of the system is now given by 

2 3 

tfcond -EE {h 9o |l>i0'l &t + \ m |0>i<2| + h.c. 

i=l j=2 

-i/sr 3 |3>i<3| . (10) 

Proceeding as above leads to the effective Hamiltonian 

H cS = ^=m[- |001) (0a\ + |010) (0a\ + h.c. ] (11) 

and if T = 2n/Sl, then the time evolution of the system 
exchanges the amplitudes of the qubit states 1 01) and 1 10) 
of the initial state. Deviations of the fidelity from unity 
are, for the same parameter regime as considered above, 
about the same size as for the proposed phase gate since 
implementation of the two gates is very similar. The 
same applies to the improvement of the gate precision 
that can be achieved with the help of dissipation and to 
the gate success rates. 

Summarising, we have shown that dissipation can be 
used to construct fast, simple and precise gates for quan- 



tum computing. As an example we discussed the imple- 
mentation of a two-qubit phase gate and SWAP opera- 
tion with cold trapped ions. Here the parameter regime 
has been chosen such that auxiliary decay channels, re- 
sulting from continuous cooling of the ions, stabilise the 
desired adiabatic time evolution of the system and im- 
prove the fidelity of the gate operation significantly. The 
more general regime where strong ground state cooling 
allows the implementation of arbitrary single laser pulse 
gates will be discussed elsewhere [13]. 
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